On the Field Theoretical Approach to the Anomalous Scaling in Turbulence 
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Anomalous scaling problem in the stochastic Navier-Stokes equation is treated in the framework of 
the field theoretical approach, successfully applied earlier to the Kraichnan rapid advection model. 
Two cases of the space dimensions d = 2 and d ^ oo, which allow essential simplification of the 
calculations, are analyzed. The presence of infinite set of the Galilean invariant composite operators 
with negative critical dimensions in the model discussed has been proved. It allows, as well as 
for the Kraichnan model, to justify the anomalous scaling of the structure functions. The explicit 
expression for the junior operator of this set, related to the square of energy dissipation operator, 
has been found in the first order of the e expansion. Its critical dimension is strongly negative in 
the two dimensional case and vanishes while d — > oo. 
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It is well known that the Kolmogorov-Obuchov (KO) 
theory fails to describe certain features of fuUy developed 
turbulence. In particular, a good example is that of the 
single time structure functions of velocity 



Snir) = {[vLir) - VLiO)]") , «L 



(vr) 



(1) 



From the viewpoint of the KO theory, their behavior in 
the inertial range depends only on the mean energy dis- 
sipation rate e 



Sn{r) = C„(£-r)"/3, 



(2) 



which leads to linearity of the corresponding exponent 
in n. However, this is not consistent with experimental 
results. This fact, often referred to as anomalous scaling, 
is being actively discussed now. The problem is usually 
treated in the framework of phenomenological models M , 
where it is related to large fluctuations of the energy dissi- 
pation rate. The deviations from the Kolmogorov scaling 
are described by additional (anomalous) exponents g„ as 
follows 



5„(r)~C„(£-r)"/3(r/i)«" 



(3) 



where L is the integral scale of the flow. For n > 3 the 
experiment shows that (/„ is notably less then zero. 

The only statistical model allowing the calculation of 
anomalous exponents is the Kraichnan rapid advection 
model, introduced by Obukhov |^ and Kraichnan [Q, 
which attracts considerable interest being a simple model 
for anomalous scaling investigation [Q. Recently it was 
treated within the held theoretical approach, based on 
the analysis of critical dimensions of so called compos- 
ite operators — products of primary field and/or their 
derivative at a single point [pj. This approach provides 
one of the most efficient ways to determine the univer- 
sal scaling quantities, by ignoring all inessential details, 
which is a great advantage in a complex problem investi- 
gation. Particularly, it became possible to calculate the 
anomalous exponents for the Kraichnan model in the sec- 
ond order of the e expansion (the flrst order result in this 



case is consistent with that obtained in [El). Another 
convenience of the approach in question is its universal- 
ity: it does not exploit the peculiarity of the Kraichnan 
model and thus can be applied to a wide variety of the 
statistical models. 

In the present letter this approach is developed for the 
stochastic Navier-Stokes (SNS) equation in two cases of 
space dimensions d = 2 and d ^ oo, which allow essential 
simplification of calculations. The obtained results are 
found to be analogous to that for the Kraichnan model. 
However, due to considerable technical complexity it is 
only qualitative analysis that can be performed. 

The starting point in the scaling investigation is the 
dimensional analysis. We will distinguish momentum d^ 
and frequency d'^ dimensions, so the dimension of an ar- 
bitrary quantity F is represented as [F] = (length) ^ x 
(time)" " . It is obvious that d^ — —d^ = 1 and 
d'^ = —d^ — 1. For the viscosity v one obtains dj; = —2, 
d'^ — 1. The fully developed turbulence has two char- 
acteristic lengths: the integral scale L, determined by 
the geometry of the flow, and the dissipative length 
I = e~^l^v^l^ . Thus, the natural representation for an 
arbitrary function S(r) is 



S{r) 



-dsp 



(I'l) 



(4) 



d% 



2d'g is so called canonical dimension of 



Here ds 

S, which plays an important role in the field theoretical 
approach (as well as conventional dimension for the static 
problems). Particularly, for the structure functions (|l|) 
d% — n, dg = —n and consequently ds„ — n, so the 



s„ 
Eq. 



takes the form 



5„(r) = z."r-"^„(^,^). 



(5) 



However, this representation is not fruitful for the fully 
developed turbulence, because of the divergences of the 
functions Fn while the Reynolds number goes to infinity. 
More precisely, one is interested in the structure func- 
tions behavior in the inertial region L 3> r 3> Z, which 
corresponds to the limit r/l — * cxd, r/L — > in Eq. (H). 



It should be noted that there are two different problems: 
the asymptotical behavior of the functions Fn while the 
first argument goes to infinity, and while the second one 
goes to zero. 

The common way for treating the first problem is to 



isolate the factor (r/Z) ^" in Eq. (||) 



Sn{r) = v^r-^ (y) 
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If the choice of the exponents 7„ can provide the con- 
dition Fn{oo,r/L) ~ const / 0, then infrared scaling of 
structure functions Sn (r) with fixed parameter r/L 

5„(r)=A^"5„(rA) (7) 

will be described by their critical dimensions A„ = ds„ + 

In- 

In the framework of the KO theory, functions 
Fn{r/L) = Fn{oo,r/L) are assumed to have nonzero 
limit while r/L —f and the exponents 7„ are determined 
from the condition of the disappearing of v dependence 
in the asymptotical behavior of the structure functions 
in the inertial region (this leads to the Eq. (pi)). 

On the other hand, critical dimensions of different 
quantities can be calculated directly from the correspond- 
ing stochastic equation by means of the renormalization 
group (RG) method. For the structure functions (|l|) the 
exact nonperturbative result can be obtained [|6|, which 
gives Kolmogorov's values A„ = — n/3. However, the 
RG method does not require the finiteness of the func- 
tions Fnir/L) in Eq. (||) (the only condition is I <S^ L). 
It is obvious that the anomalous scaling appears in the 
case of divergencies of the functions F„ while r/L -^ 0. 

Within the field theoretical approach, the asymptoti- 
cal behavior of the functions Fn while r/L — > can be 
analysed with the help of the operator product expan- 
sion (OPE), which allows to express the dependence dis- 
cussed in terms of critical dimensions of composite opera- 
tors. The detailed account for the RG and OPE methods 
applied to SNS equation can be found in |q|, so here we 
provide only necessary information. In general, OPE of 
the discussed above functions Fn has the form 



F„cx^CW(r/L)^- 



(8) 



where the coefficients C„ are regular functions of r/L 
and thus can be treated as constants in the limit r/L -^ 
0. Summation in Eq. (S) is taken over all composite op- 
erators $fc, entering the OPE of 5„ (in general, these 
are all possible operators allowed by the symmetry) , and 
A$^, are their critical dimensions. The latter, unlike the 
critical dimensions A„ of the structure functions, cannot 
be determined from the simple dimensionality consider- 
ations. In this case the application of the RG method 
is required, which allows to calculate nontrivial critical 



dimensions A$j. as well as coefBcients C„ in the form 
of e expansions. 

One can see from Eq. (g|) that the main contribution 
to Fn in the asymptotic region r/L ^ is given by terms 
with negative A^, i. e. by composite operators with neg- 
ative critical dimensions. These operators were called 
"dangerous" in , because the presence of such opera- 
tors in the expansion (g) leads immediately to the diver- 
gences of the functions F„ while r/L -^ and, hence, 
to the anomalous scaling. As it was pointed out in [||, 
only Galilean invariant composite operators enter into 
the OPE of invariant objects, as the structure functions 
are. But it was not until now that invariant dangerous 
operators have been found in the SNS model. It is only 
the energy dissipation operator 



^d = -^J^idiifk + dkipif , (p = v 



(v), 



(9) 



that was known to have zero critical dimension [|6| . 

One of the main difficulties of the search of dangerous 
operators in the SNS model originates from the mixing 
of operators in renormalization. The matter is that an 
arbitrary ultraviolet finite renormalized operator is a lin- 
ear combination of the primary one and those mixing in 
renormalization. This results in splitting of the whole set 
of operators into finite subsets — "families" , closed un- 
der renormalization and containing operators of the same 
canonical dimensions and symmetry properties |q| . As a 
consequence, the renormalization of a composite opera- 
tor <I>i is characterized by a renormalization matrix Zij, 
rather than by a single constant as in the case of ordinary 
quantities 



$, 



E^^.'J'f- 



(10) 



Here the sum is taken over the full set of independent op- 
erators from the corresponding operator family. In fact, 
definite critical dimensions may be assigned to the so 
called basis operators only, which diagonalize the matrix 
Zij, and are determined by its eigenvectors. 

Thus, the RG analysis of an operator family re- 
quires calculating of the full renormalization matrix from 
Eq. (IG). This provides possibility to determine all basis 
operators and their critical dimensions. But the number 
of independent operators in the SNS model grows rapidly 
with the canonical dimension of the family, which makes 
the analysis of senior families extremely laborious. 

Particularly, we will be interested in the scalar Galilean 
invariant operators as entering the OPE of the structure 
functions. Canonical dimensions of the operators in ques- 
tion are positive even integers. Up to the present moment 
the families with (i$ = 2,4,6 have been examined [^,||. 
The mentioned above energy dissipation operator ^^ be- 
longs to the family with d$ = 4, and has the minimal 
critical dimension A$^ = among the families discussed. 



This fact encourages us to consider the famihes, con- 
taining powers of $d, to be the most hkely evidence of 
the presence of dangerous operators. The family with 
d$ = 8, which contains the square of <i>d, was examined 
in pOJ. However, the analysis of the whole family was not 
performed due to technical difficulties (the family con- 
tains 12 independent operators). It was found out only 
that <i>^ cannot be a basis operator because of its mixing 
with the others: the fact which breaks the equality 



A[$3]=nA[$rf], 



(11) 



proposed in |lll| . 

We succeeded in the analysis of the discussed above 
family with canonical dimension 8 in two cases of space 
dimensions d = 2 and d — > oo, where the essential simpli- 
fication of the problem occurs. Let us consider the first 
case in more detail. 

In this case the SNS equation for the velocity pulsation 
field (fi can be reduced to the scalar equation 



dtAcj) = l^A'^cl) - did„i{£nmdn(t)di> 



f 



(12) 



for the stream function 



'jiOjC, 



defined by the relation ipi — 
2 13|. Here e,j is the second-rank antisymmetric 
tensor. Correlator of the random force /, which imitate 
the interaction of velocity pulsations with large-scale ed- 
dies, can be taken in the usual form 

{fix,t)fix',t'))^S{t-t') 

^ exp {zk(x - x')} i?(k) , (13a) 



(27r)2 



D(k) = Dak 



4-2e 



where e > is a formal small parameter of RG expan- 
sion. Its "physical" value e = 2 corresponds to the case 
of energy pumping. Correct RG approach for this model, 
developed in O-Q gives the same results, as the stan- 
dard SNS model. 

Scalar formalism has at least two advantages: simpler 
calculation of each diagram and less number of indepen- 
dent operators. Let us see how the family of operators 
with canonical dimension d$ = 8 can be treated within 
the scalar approach of the model (|lj) , (|l^) . To analyse 
an operator family one has to found the full set of in- 
dependent operators. Operators should be constructed 
from the derivatives 9^, dt of field convoluted with 
tensors £jj, Sij and must necessarily be true scalars (an 
essential demand as Sij is a pseudo-tensor and is a 
pseudo-scalar) and Galilean invariant. It can be shown 
that the consideration may be restricted to the case of op- 
erators irreducible to a full derivative of a junior one M. 
Also, the specific symmetries of two-dimensional space 
must be taken into account, leading to identities of the 
type SijEki = SikSji - SiiSjk- 



As the result one obtains 7 independent operators: 
$1 = (A0)4, 

^3 = didj4> d.jdk4> dkdicj) didi4>, 



$4 = £ii Adidjcj) djdkct) dkdicj) 
$7 = Vtdidj<j)\7td^dj(f>. 



(14) 



Here Vt = 9t + ejidjcjidi is the substantial derivative. 
The square of energy dissipation operator is expressed in 
terms of operators (|lj) as follows 



<I>^ = -$1 + 6<i>2 -4<I>3 



(15) 



We have calculated full 7x7 rcnormalization matrix 
Zij for the set of operators (|lj) in the first order of the e 
expansion within the minimal subtraction scheme. As it 
was claimed in [Q , <I>^ is not a basis operator and does 
not have any certain critical dimension. But a dangerous 
operator still enters into the family in question. It is 



*'.^ 



$1 + 6.03$2 - 4.04$3 - 0.58$4 
f 0.06$5 + 0.10$6-0.02$7 



with critical dimension 



A 



*'. 



8 - 5.92e. 



(16) 



(17) 



At the "physical" point e = 2 it becomes A$' = —3.84. 
It is interesting to notice that the square of energy dis- 
sipation operator (nS) makes the major contribution to 
the dangerous operator <I>'^, while the other operators of 



the family (14) enter with small numerical coefficients. 



(13b) 'pj^g same is the case with the critical dimension of $ 



the major contribution 8 — 6e is that of diagonal clement 
of the renormalization matrix, corresponding to the <I>|, 
and the admixtion correction comprises only 0.08e. So 
the operator (|l^) proves to be dangerous due to the 
square of the energy dissipation operator. However, it 
must be outlined that the discussed critical dimension 
A$' considerably deviates from 2A$^ = 8 — 4e given 
by(0). 

It is the second studied case of infinite dimensional tur- 
bulence, where the exclusive role of the energy dissipation 
operator becomes even more obvious. We have shown 
that this is the case when the renormalization matrix is 
block-triangular, so that the <I>^ operator is the basis one. 
Its critical dimension 



A<B2 = 8 - 4e 



(18) 



satisfies the Eq. ( pi] ) and equals to zero in the "physi- 
cal" point as for the energy dissipation operator itself. 
Moreover, we show within one-loop approximation that 
in the case discussed the Eq. ww proves to be valid for 
all powers of <i>d p^] . 



Let us show now that the existence of one dangerous 
operator in the expansion (H) means that there is an in- 
finite set of dangerous operators. Indeed, by combining 
inequaUty 



l^nl'/" < \Sn+l 



|l/(n+l) 



(19) 



known from the probabihty theory with the definition of 
anomalous exponents (0) , one obtains inequahty for even 
exponents <?„ 



On 

n 



< 



Qn+2 



2m. 



(20) 



Since qn are negative (which foUows from both the exper- 
imental data and the presence of a dangerous operator 
in the OPE of structure functions) the inequality (Eo) 
implies that the absolute value of qn grows more rapidly 
than n. Taking into account the expansion m) we see 
that only infinite set of dangerous operators meet the 
last condition. 

Thus, the results, obtained for the SNS model, proved 
to be quite similar to those for the Kraichnan model [Q . 
In both cases the anomalous scaling of the structure func- 
tions can be considered as a consequence of the infinite 
set of dangerous operators in the corresponding OPE, 
related to the powers of the energy dissipation opera- 
tor. Within the used approach, both models exhibit 
the particularity of the infinite-dimensional case, were 
the critical dimensions of all known dangerous opera- 
tors vanish (at least within the considered approxima- 
tions). In the Kraichnan model this is in agreement with 
the disappearing of anomalous scaling and intermittency 
while d ^ oo |16). Our results make us expect that the 
analogous behavior takes place in the SNS model also, 
which was first discussed in pTJ, where the possible use 
of 1/d expansion for real turbulence was pointed out. 

However, the Kraichnan model contains at least two 
principal simplifications. Firstly it is free of problems 
connected with the mixing of operators: the critical di- 
mension of any power of the dissipation operator is de- 
termined by the corresponding diagonal element of the 
renormalization matrix (due to its triangularity). This 
allows to calculate the critical dimensions of all opera- 
tors discussed in e expansion |g| , which seems impossible 
for the SNS model. Secondly, the expansion (|8[) for any 
certain structure function Sn is finite in the Kraichnan 
model, which eliminates the necessity of its summation. 
The anomalous exponent (/„ in this case simply equals to 
the minimal Aj, in Eq. (|^). However there is no reason 
to expect this in the SNS model, so the additional prob- 
lem of the summation of the expansion (g) arises. As a 
consequence, the asymptotical behavior of the structure 
functions can be more complicated than powerlike. 

It becomes clear now that the developed here approach 
gives qualitative account for anomalous scaling in the 
SNS model, though, due to extreme difficulty of the prob- 



lem it requires some additional ideas for further investi- 
gations. In this context the simple result for the infinite- 
dimensional case is of particular interest. It implies the 
actuality of the detailed study of this case as the starting 
point of 1/d expansion. 
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